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In the framework of seesaw mechanism embedded in the constrained Minimal Supersymmetric 
Standard Model (cMSSM), phases of neutrino Yukawa coupling, /i-term and ^4-terms can all con- 
tribute to the Electric Dipole Moment (EDM) of the electron. We discuss and classify the situations 
for which by combined analysis of the upcoming results on d e , dn s and do discriminating between 
these sources will be possible. 



I. INTRODUCTION 

Recent atmospheric and solar neutrino data as 
well as the results of KamLAND 0, K2K and MI- 
NOS establish nonzero mass for neutrinos. On the 
other hand, kinematical studies fa] and cosmological 
observations show that neutrino masses cannot ex- 
ceed a few eV. Nonzero neutrino masses cannot be 
accommodated within the Standard Model (SM) . Sev- 
eral models have been suggested to attribute tiny yet 
nonzero masses for neutrinos among which the see- 
saw mechanism [?| is arguably the most popular one. 
The seesaw mechanism adds three SM singlet right- 
handed neutrinos with very heavy Majorana masses to 
the model. Up to now, we have no clue how large the 
masses of the right-handed neutrinos are. The scale 
can lay anywhere between TeV up to 10 14 GeV. If the 
scale is higher than 1 TeV, the right-handed neutrinos 
cannot be produced by accelerator technology; thus, 
we can only learn about the seesaw parameters indi- 
rectly through their effects on low energy parameters 
0, such as light neutrino mass matrix and slepton 
mass matrix in the context of supersymmetric seesaw 
0. In this line, any low energy observable which is 
sensitive to the seesaw parameters deserves special at- 
tention. 

The 3x3 neutrino Yukawa matrix introduces six 
new sources of CP-violation which, in the context of 
supersymmetric standard model, can induce signifi- 
cantly large contributions to the electric dipole mo- 
ments of the electron lllj and ot her p articles. So 
far no finite EDM has been observed [12|: 

|4| < 1.4 x 10~ 27 e cm. (1) 

However, there are ongoing experiments [l3| as well as 
proposals [Llj to improve the present bound by several 
orders of magnitude. In view of these experiments, it 
has been suggested to use the EDMs to extract infor- 
mation on the seesaw parameters [15J. However, for 
deriving any information from d e we must be aware of 
other sources of CP-violation that can give a signifi- 
cant contribution to d e . 

In the cMSSM, there are two extra sources of CP- 
violation relevant for the EDMs of leptons: the phases 



of the fj, parameter and the universal trilinear coupling 
ao. The phases of ao and /i induce d e ~ (m e /md)dd ~ 
(m e /m u )d u ~ e(m e /m d )dd ~ e(m e /m u )d u , where d u 
and dd, respectively, are the Chromo EDMs (CEDM) 
of up and down quarks which give rise to EDMs of 
hadrons and nuclei such as mercury and deuteron. 
However, the quark EDMs and CEDMs induced by 
the phases of Y v are too small to be detectable in near 
future. Therefore, if complex Y v is the only source 
of CP-violation, we expect the Deuteron EDM (do) 
to be too small to be detectable in the near future 
{do is measured with ionized deuteron which is de- 
pleted from electrons). Based on this difference, it 
has been suggested in |16| to combine the informa- 
tion on do and dn g with d e to disentangle the source 
of CP-violation. It is also discussed how much the 
present bounds have to be improved in order to be 
able to make such a discrimination. In present letter, 
we review the results obtained in Ha. 



II. THE MODEL 

The seesaw mechanism embedded in the MSSM is 
described by the superpotential 

W = Y-E.Hd-L^N^.L^M^N^Hd-H, 

(2) 

The quark Yukawa couplings, not shown here, are the 
same as in the MSSM. Here, i,j are generation in- 
dices, Lj consist of lepton doublets (^l^Jl), an d Ei 
contain left-handed anti-leptons if L , The superfields 
Ni contain anti right-handed neutrinos. Without loss 
of generality, one can rotate and rephase the fields 
to make Yukawa couplings of charged leptons (Ye) as 
well as the mass matrix of the right-handed neutrinos 
(Mij) real diagonal. In what follows, we will use this 
basis. 

In general, the soft supersymmetry breaking terms 
(the mass-squared matrices and trilinear couplings 
of the sfermions) can possess flavor-changing entries 
which facilitate a number of flavor-changing neutral 
current processes in the hadron and lepton sectors. 
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The existing experimental data thus put stringent 
bounds on flavor-changing entries of the soft terms. 
For instance, flavor-changing entries of the soft terms 
in the lepton sector can result in sizeable fj, — > &y, 
t — ► B*y and r — > [ij. This motivates us to go to 
the constrained MSSM framework where soft terms 
of a given type unify at a scale close to the scale of 
gauge coupling unification. In other words, at the 
GUT scale, we take 

Csoft = -mg(Z i t L< + E\ Ei + N}Ni + HlH d + HlH u ) 

- ~m 1/2 (BB + WW + gg + Tl.c.) 

- (^B H H d -H u + R.c.) 

- (4 j EiH d ■ Lj - A^NiH u ■ Lj + H.c.) 

- (iB tf M i JV*JV i + H.c). (3) 

Here Ag = aoYf and A v = a§Y v , where ao in general 
can be complex and a source of CP-violation. The last 
term is the lepton number violating neutrino bilinear 
soft term which is called the neutrino B-term. As has 
been first shown in 0] , the phase of the neutrino B- 
term can induce a contribution to d e . In this letter, 
for simplicity, we will set B v = 0. By rephasing the 
Higgs fields we can make Bh real; however, the phase 
of n will in general remain nonzero. In this letter, to 
calculate the effects of the phases of Y v , \x and ao on 
the EDMs and CEDMs, we use the results of [TlllTl 
0. 

III. BOUNDS ON THE EDMS 

In this section, we first review the current bounds 
on do-, dug and d n and the prospects for improving 
them. We then review how we can write them in terms 
of Im(/i) and Im(ao). 

• Neutron EDM, d n : The present bound on d n 
[His 

d n < 3.0 x l(T 26 e cm at 90% C.L. (4) 

This bound can be improved considerably by 
SNS [2(| which will be able to probe d n down to 
10~ 28 e cm. 

• Mercury EDM, d# g : The present bound on 
d Hg is 

\d Hg \ < 2.1 x 10~ 28 e cm (5) 
which can be improved by a factor of 2 pi) . 

• Deuteron EDM, do- The present bound on 
do is too weak to constrain the CP-violating 
phases; however, there are proposals |22( to 
probe do down to 

(1 - 3) x 10~ 27 e cm. (6) 



Different sources of CP-violation affect the EDMs 
listed above differently. As a result, in principle, 
by combining information on these observables we 
can discriminate between different sources of CP- 
violation. It is rather straightforward to calculate 
the EDMs and CEDMs of the elementary particles in 
terms of the phases of ao and fj, 0,0,01; however, 
writing d„, dn g and do in terms of the EDMs and 
CEDMs of their constituents is more difficult and a 
subject of debate among the experts. Let us consider 
them one by one. 

• d n (d q ,d q ): 

Despite the rich literature on d„ in terms of the 
quark EDMs and CEDMs, the results are quite 
model dependent. For example, the SU(3) chiral 
model 23] and QCD sum rules 0] predict dif- 
ferent contributions from d u and dd to d n . Con- 
sidering these discrepancies in the literature, we 
do not use bounds on d n in our analysis. 

• d Hg {d qi d q ): 

There is an extensive literature on du g 0] . Fol- 
lowing Ref. p^ ] , we will interpret the bound on 
d Hg as 

\d d -d u \ < 2 x 10" 26 cm. (7) 

• d D {d q ,d q ): 

Searches for do can serve as an ideal probe for 
the existence of sources of CP-violation other 
than complex Y v because i) there is a g ood 
prospect of improving the bound on do |22j : ii) 
an ionized deuteron does not contain any elec- 
trons and hence we expect only a negligible and 
undetectable contribution from Y v to do- 

To calculate do in terms of quark EDMs and 
CEDMs, two techniques have been sugg ested in 
the literature: i) QCD sum rules [27| and ii) 
SU(3) chiral theory 0]. Within the error bars, 
the two models agree on the contribution from 
dd~d u which is the dominant one. However, the 
results of the two models on the sub-dominant 
contributions are not compatible. Apart from 
this discrepancy, there is a large uncertainty in 
the contribution of the dominant term: 

d D {d q ,d q ) ~ -e(d u - d d ) . (8) 

In this paper we take "the best fit" for our anal- 
ysis. 

IV. NUMERICAL ANALYSIS 

In this section, we first describe how we produce the 
random seesaw parameters compatible with the data. 
We then describe figs. (Ill 111 and, in the end, discuss 
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what can be inferred from the future data considering 
different possible situations one by one. 

In figures (111 111 , the scatter points marked with "+" 
represent d e resulting from complex Y u . To extract 
random Y v and Mjv compatible with data, we have 
followed the recipe described in and solved the 
following two equations 

Yj— Y v (v 2 sin 2 0)/2 = U • $ • M® iag • $ • U T (9) 



and 



v+T M GUT V 



a d 
6 
d* c 



(10) 



where v — 247 GeV, M is the mass matrix of the right- 
handed neutrinos, U is the mixing matrix of neutrinos 
with S13 = and $ is diag[l, e 1 ^ 1 , e 1 ^ 2 ] with random 
values of <j>\ and <jf>2 in the range (0, 2w). Finally, 
diag[mi, \J m\ + Am^, y^r 



M m ag 



where mi picks up random values between and 0.5 
eV in a linear scale. The upper limit on mi is what has 
been found in by taking the dark energy equation 
of state a free (but constant) parameter. 

In order to satisfy the strong bounds on Br(/i — * ej) 
[l2| and Br(r — > ^7) jU, the matrix ft,, defined in Eq. 
Q10JI. is taken to have this specific pattern with zero 
e/x and \it elements. Actually these branching ratios 
put bounds on (Am? ) eM and (Am? ) jUT rather than on 
/i e /n and /i^t-. Notice that only the dominant term of 
Am? is proportional to h. There is also a subdomi- 
nant "finite" contribution to Am? which is about 10% 
of the dominant effect and is not proportional to the 
matrix h Nonetheless, for extracting the see- 

saw parameters, 20% accuracy is enough and we can 
neglect the subdominant part and take Am? propor- 
tional to the matrix h. In Eq. (|1C)|> . a, 6, c are real 
numbers which take random values between and 5. 
On the other hand, \d\ takes random values between 
and the upper bound from Br(r — > ej) ,32]. To 
calculate the upper bound on \d\, we have used the 
formulae derived in Ref. j3^|. The phase of d takes 
random values between and 2ir. 

To perform this analysis we have taken various val- 
ues of tan (3 and ao and calculated the spectrum of the 
supersymmetric parameters along the m^- mp strips 
parameterized in Ref. [34j]. Notice that Ref. [34[ has 
already removed the parameter range for which color 
or charge condensation takes place. 

In the figures, we have also drawn the present bound 
on d e [T2]] as well as the limits which can be probed in 
the future. The present bound is shown by a dashed 
dark blue line and lies several orders of magnitude 
above the d e from phases of Y v . After five years of 
data-taking, the Yale group can probe d e down to 
10~ 31 e cm 0| which is shown with a dot-dashed 
cyan line in the figures. As demonstrated in the fig- 
ures, only for large values of ao the effect of complex 



Y v on d e can be probed by the Yale group and for most 
of the parameter space the effect remains beyond the 
reach of this experiment. 

There are proposals [l4j to use solid state tech- 
niques to probe d e down to 10 -35 e cm (shown with 
dot-dashed yellow line in the figure). In this 
it can be deduced from the figure, we will have a great 
chance of being sensitive to the effects of the phases of 
Y v on d e . However, unfortunately, the feasibility and 
time scale of the solid state technique is still uncertain. 

In figs, ^and d e resulting from Im[/i] is also de- 
picted. The red solid lines in these figures show d e 
from Im[/i] assuming that the corresponding dn g sat- 
urates the present bound |2l| . As is well-known, there 
are uncertainties both in the value of rrid |l2j and in 
the interpretation of dn g in terms of more fundamen- 
tal parameters d M , da and d s . To draw this curve we 
have assumed m,j = 5 MeV and d u — dd < 2 x 10 -26 
e cm . As shown in the figure this bound is weaker 
than even the present direct bound on d e . The purple 
dotted lines in figs. 0}, represent d e induced by 
values of Im[/i] that give rise to d u — dd — 2 x 10 -28 
cm (corresponding to dp = 10 -27 e cm and do = 
5e(dd — d u j). Notice that these curves lie well below 
the direct bound on d e but the Yale group will be able 
to probe even smaller values of d e . Similarly in figs. 
1|2I3[1 . d e resulting from Im[a ] is depicted. 

In the figures, the bounds from du g and do appear 
almost as horizontal lines. This results from the fact 
that for the tuq — mi/ 2 strips that we analyze, mo 
is almost proportional to m 1( / 2 . Using dimensional 
analysis we can write 



fci 



Im[/x] or Im[ao] 

m ?/2 



h 



Im[/i] or Im[ao] 

m \/2 



where k\ and ki are given by the relevant fermion 
masses and are independent of m^ 2 . As a result, for 
a given value of d u — dd, Im[/i] (or Im[ao]) itself is 
proportional to m 3 ^ 2 so d e will not vary with m^. 

In the following, we will discuss what can be inferred 
about the sources of CP-violation from d e and do if 
their values (or the bounds on them) turn out to be in 
certain ranges. According to Fig.^ for a o = 0, any 
signal found by the Yale group implies that there are 
sources of CP-violation other than the phases of the 
Yukawa couplings. However, for larger values of ao, 
the effect of Y v on the EDMs can be observed by the 
Yale group within five years. According to Figs. 121 II) . 

for ao ~ 1000 GeV EDMs originating from complex 
Y v can be large enough to be observed by the Yale 
group. Therefore, if after five years the Yale group 
reports a null result, we can derive bounds on certain 
combinations of seesaw parameters and ao. At least 
it will be possible to discriminate between low and 
high ao values. However, if the Yale group finds that 
10~ 31 e cm < dp < 10 -29 e cm we will not be able to 
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FIG. 1: Electron EDM for a = 0, tan/3 = 10 and 
sgn(^i) = +. The scatter points represent |d e | induced 
by random complex Y u compatible with the data. The 
blue dashed line is the present bound on d e [T3 | and dot- 
dashed lines show the values of d e that can be probed in 
the future Q 
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FIG. 2: Electric dipole moment of the electron for an = 
1000 GeV, tan/3 = 10 and sgn(^) = +. To draw the red 
solid and purple dotted lines, we have assumed that Im[an] 
is the only source of CP-violation and have taken dd — d u 
equal to 2 x 10 -26 cm and 2 x 10 -28 cm, respectively to 
derive Im[ao]. To produce the scatter points, we have as- 
sumed that complex Y v is the only source of CP- violation 
and have randomly produced Y v compatible with the data. 
The blue dashed line is the present bound on d e [12] and 
dot-dashed lines show the values of d e that can be probed 
in the future [TUT^I 



FIG. 3: The same as fig. EJfor a = 1000 GeV, tan = 20 
and sgn(^i) = +. 



determine whether d e originates from complex Y v or 
from more familiar sources such as complex ao or ji. 
To be able to make such a distinction, values of 
down to 1CT 28 - 1CT 29 e cm have to be probed which, 
at the moment, does not seem to be achievable. 

If future searches for do find dp > 10" 27 e cm but 
the Yale group finds d e < 2 x 10~ 29 e cm (this can be 
tested within only 3 years of data taking by the Yale 
group ^3) j we might conclude that the source of CP- 
violation is something other than pure Im[/i] or Im[ao]; 
e.g., the QCD 9-term which can give a significant con- 
tribution to dp but only a negligible contribution to 
d e . Another possibility is that there is a cancelation 
between the contributions of Im[/i] and Im[ao] to d e 
|35j]. The information on d n would then help us to 
resolve this ambiguity provided that the theoretical 
uncertainties in calculation of d n as well as are 
sufficiently reduced. 

On the other hand, if the Yale group detects d e > 
2 x 10~ 29 e cm, we will expect that d D > 10~ 27 
e cm which will be a strong motivation for build- 
ing a deuteron storage ring and searching for djj. If 
such a detector finds a null result, within this frame- 
work the explanation will be quite non-trivial requir- 
ing some fine-tuned cancelation between different con- 
tributions. 

According to these figures, in the foreseeable future, 
we will not be able to extract any information on the 
seesaw parameters from EDMs, because even if we 
develop techniques to probe d e as small as 10~ 35 e cm, 
we will not be able to subtract (or dismiss) the effect 
coming from Im[/i] and Im[a ] unless we are able to 
probe d u — dd at least 5 orders of magnitude below its 
present bound which seems impractical. Remember 
that this is under the optimistic assumptions that the 
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FIG. 4: Electric dipole moment of the electron for an = 
2000 GeV, tan/3 = 20 and sgn(^) = +. To draw the red 
solid and purple dotted lines, we have assumed that Im[/j] 
is the only source of CP-violation and have taken dd — d u 
equal to 2 x 10 -26 cm and 2 x 10 -28 cm, respectively to 
derive Im[/j]. To produce the scatter points, we have as- 
sumed that complex Y v is the only source of CP- violation 
and have randomly produced Y u compatible with the data. 
The blue dashed line is the present bound on and 
dot-dashed lines show the values of d e that can be probed 
in the future (HQ 

mass of the lightest neutrino, mi, and Br(r — > ej) are 
close to their upper bounds and there is no cancelation 
between different contributions to the EDMs. 

If, in the future, we realize that mi and Br(r — > 
ej) are indeed close to the present upper bounds on 
them and ao = (oo = 1000 GeV) but find d e < 
10~ 35 e cm (d e < 10~ 34 e cm ), we will be able to 
draw bounds on the phases of Y v which along with 
the information on the Dirac and Majorana phases of 
the neutrino mass matrix and the CP-violating phase 
of the left-handed slepton mass matrix may have some 
implication for leptogenesis. This is however quite an 
unlikely situation. 

V. SUMMARY 

In this work we have studied EDMs of particles in 
the context of supersymmetric seesaw mechanism. In 



figs. Ill II the values of d e corresponding to differ- 
ent random complex Y v textures are represented by 
"4-". For small values of tan/3 (tan/3 < 10) and ao 
(ao < 1000 GeV), d e induced by Yw is beyond the 
reach of the ongoing experiments 13]. Such values 
of d e can however be probed by the proposed solid 
state based experiments . For larger values of tan /3 
and/or a , the Yale group may be able to detect the 
effects of complex Y v on d e . As demonstrated in Figs. 
01 and H for tan [3 = 20 and a = 1000 - 2000 GeV, a 
large fraction of parameter space yields d e detectable 
by the Yale group. However, even in this case we will 
not be able to extract information on the seesaw pa- 
rameters from d e because the source of CP-violation 
might be ao and/or fj, rather than Y v . If the future 
searches for do |23 find out that do > 10~ 27 e cm 
then we will conclude that there is a source of CP- 
violation other than complex Y v . However, to prove 
that the dominant contribution to d e detected by the 
Yale group comes from complex Y v — hence to be able 
to extract information on the seesaw parameters from 
it- we should show that do < 10~ 28 — 10~ 29 e cm 
which is beyond the reach of even the current pro- 
posals. Notice that for the purpose of discriminat- 
ing between complex Y v and ao/fi as sources of CP- 
violation, searching for dn g is not very helpful because 
mercury atom contains electron and hence dHg ob- 
tains a contribution from complex Y v . That is while 
ionized deuteron used for measuring du does not con- 
tain any electron and the contribution of complex Y v 
to it is negligible. To obtain information from d n , the 
theoretical uncertainties first have to be resolved. 
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